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Abstract 

We have carried out density functional theory based calculations of the size dependent formation 

energy and geometry of Pt islands on Ru(OOOl), to model Pt-Ru nanocatalysts which have been 

recently proposed as fuel cell anode. The Pt islands are found to prefer two-dimensional structures. 

Furthermore, a monotonic decrease in the formation energy per Pt atom suggests a propensity of 

^ I Pt atoms to wet the Ru(OOOl) surface. Calculated energy barriers for the diffusion of Pt monomers 

'^ ■ and dimers on the facets and through the edges of a superstructure modeling a Ru nanoparticle 

O . indicate that these edges help reduce considerably the diffusion rates across them such that the Pt 

atoms prefer to remain in the facet on which they were adsorbed originally and form 2D islands. 

(N 
> 

Q> \ PACS numbers: 

T-H ' 

in 

i> 

o 

00 

O 
> 

X 



I. INTRODUCTION 

Direct methanol fuel cells (DMFC) are considered a promising means for energy conver- 
sion in "hydrogen-based economy" because they work at low temperatures (~ 350-400 K) 
and use liquid methanol as fuel, which is easy to deliver and store. In DMFC, the same 
anode is used as a catalyst for both methanol reforming and for the oxidation of hydrogen 
obtained from the reforming. Although the carbon monoxide released in the course of this 
reaction is expected to be oxidized by hydroxyl radicals obtained from admixed water, ex- 
periments find that it still severely poisons the commonly used Pt anode by blocking the 
reactive Pt sites and, thereby, reducing the rate of hydrogen electro-oxidation. Similarly, in 
proton exchange fuel cells operating with pre-reformed gas, the anode is poisoned by carbon 
monoxide molecules, inevitably present in hydrogen obtained from hydrocarbons. 

It is known- that PtRu alloys are more tolerant to CO poisoning than pure Pt, though 
their tolerance is still unsatisfactory. The high content of expensive platinum in these alloys 
also speaks against their choice as cost effective catalysts. It is thus encouraging to note 
that nanoclusters of Ru exposed to spontaneous Pt deposition are much more tolerant to 
CO than commercial PtRu catalysts,-*^ particularly since the content of Pt on these novel 
materials is significantly smaller than that in Ru-Pt alloys. For example, Brankovic et al.—'^ 
deposited Pt on ~2.5 nm size Ru nanoparticles and found that the 1:20 ratio (PtRu2o), 
which corresponds to ~ 0.1 monolayer (ML) coverage, surpass substantially the catalytic 
performance of PtRu and Pt2Ru3 in the presence of CO. Assuming that Pt atoms form 
two dimensional (2D) islands on the facets of the Ru nanoparticles,^ the CO tolerance was 
attributed^ and explained^ on the basis of the spillover effect, which refers to the tendency 
of CO to leave the Pt island and diffuse towards Ru. 

As noted above, Pt coverage is critical for the catalytic properties of the Ru nanopar- 
ticles.-)^ Island size effects - tuned by Pt coverage - have also been reported for methanol 
electro-oxidation.- Earlier calculations have indicated that the CO adsorption energies on 
binary Pt-Ru systems, ranging from Pt(lll) to surfaces of ordered Pt-Ru alloys to IML of 
Pt on Ru(OOOl), are reduced upon lowering Pt content.- Moreover, our recent calculation^ 
supporting the idea of CO spillover relies on the presence of small 2D Pt islands on the facets 
of Ru nanoparticles. All of the above suggest that the correlation between catalytic activity 
and coverage has implications on how Pt arranges itself on the Ru nanoparticles, as we shall 



see later. Indeed, the first step in obtaining a systematic understanding of the enhanced 
reactivity of this catalyst is the determination of the geometry and relative stability of the 
Pt islands on the nanoparticle facets. Towards this end we have adopted in this work a 
theoretical approach based on ab initio electronic structure calculations. 

Ever since the work by Stranski in 1927 about crystal growth, a large amount of effort 
has been dedicated to understanding the microscopic processes that control nucleation and 
epitaxial growth.- It is now recognized that the resulting structure when adatoms are ad- 
sorbed onto a surface is thermodynamically and kinetically controlled. The energetic barriers 
for the various processes to occur, the surface and/or post-growth annealing temperature, 
along with the deposition rate and technique, essentially determine the morphology of the 
deposited adatom structure at a given coverage. Three major types of growth modes have 
been observed in experiment.- The Frank van der Merwe mode refers to the two dimensional 
(2D) layer-by-layer growth, the Volmer- Weber growth describes three-dimensional (3D) clus- 
tering of adatoms on the bare substrate, while the Stranski-Krastanov mode consists of a 
mixture of the above two growth modes, i.e., 3D clustering occurs after one or a few com- 
plete adlayers are formed.-'^ Theoretical models of heteroepitaxial growth suggest that the 
growth mode is determined by the competition of factors such as the surface energies of 
the bare substrate and the heteroepitaxial layer, the interface free energy, and the strain 
energy introduced by the lattice mismatch of the two species.-i^ii^ For example, a mismatch 
in the respective bulk lattice constants strains the interface and may set off the 3D cluster- 
ing growth mode.- On the contrary, growth of adlayers with lower surface energy than the 
substrate may favor the 2D layer-by- layer growth mode.^ Considerations along these lines, 
however, lead to ambiguities in predicting Pt growth on Ru(OOOl). While the higher cohe- 
sive energy of Ru (relative to Pt) may^iii^ imply that the surface free energy of Ru(OOOl) 
is higher than that of Pt(lll) and point to 2D layer- by-layer growth, the stress caused by 
the Ru-Ru, Ru-Pt, and Pt-Pt bond length misfit (The bond lengths in bulk Ru and Pt are 
2.706 A and 2.775 A, respectively.—) may lead to 3D clustering. Of course, there may be 
a competition between the above two factors, leading to a critical Pt island size at which 
there is crossover between 2D and 3D growth mode or island-substrate atom exchange.— ^^^ 
One of the goals in this work is to determine whether there is indeed a critical size beyond 
which 2D Pt islands are no longer stable on Ru(OOOl). 

Turning to the case of the Ru nanoparticles described in the experiments in question,-i^ it 



is important to note that, while the experimental determination of the structure and chemical 
ordering of Pt-Ru nanoalloys represented a challenge^^ii^ii^ until the recent work of Maillard 
et al.,— the dominant features observed via surface X-ray scattering, scanning tunneling mi- 
croscopy, Fourier transform infrared spectroscopy, and high resolution transmission electron 
microscopy techniques, in 2-3 nm pure Ru nanoparticles on a carbon substrate have been 
consistent with the hexagonally close-packed (hep) Ru single crystal structure.-'i^ The Ru 
nanoparticles present well-defined facets with the lowest-surface-energy geometries, such as 
(0001), (1101), etc-ii^ In the present study, as a first step in the modeling of the Pt-decorated 
Ru nanoparticles, we focus our attention on the formation of Pt islands as a function of size 
on Ru(OOOl), one of the most predominant and stable facet orient at ions .-'^^'22 To this end, we 
carry out first principles calculations of the system's total energy to determine the geometry 
and formation energy of Pt islands, as well as that of IML of Pt on Ru(OOOl). Clearly, the 
main drawback of such calculations is that the predictions are relevant to zero temperature 
and samples relaxed in infinite time. To partially include the system dynamics, we have 
taken into account the diffusivity of Pt adatoms since, particularly at temperatures well be- 
low room temperature, the three aforementioned growth modes can be understood in terms 
of diffusion barrier differences between on-step and step-descending hopping, Schwoebel bar- 
riers.— 122>2^ It is argued that if the energy barrier at the step edge is higher (reflective barrier) 
than that of adatom diffusion on the step - positive Schwoebel barrier - there is a probability 
that adatoms will be trapped on the step terraces and 3D clustering will be favored. On 
the contrary, if the Schwoebel barrier is zero or negative (non- reflective barrier), adatoms 
that happen to lie on the step terrace are more likely to hop onto the substrate and favor 
the 2D growth.— Moreover, in heteroepitaxial growth, the mixed Stranki-Krastanov mode 
dominates if the Schwoebel barrier is reflective for homo-steps (formed by the adatoms), but 
non-reflective for hetero-steps.— 12^ 

Furthermore, since Ru nanoparticles present facets with various geometries, it is rea- 
sonable to assume that the equilibrium shape is a polyhedron and, therefore, there are 
edges delimiting the facets. Theoretical^*^ and experimental^^ studies of the shape of Pt 
nanoparticles supported on carbon, for instance, have consistently concluded that they are 
in the form of cubo-octahedra, which is a shape characteristic of fee packing. For hep 
metals, in turn, one of the proposed structures for sufficiently large nanoparticles (with 
more than 500 atoms) is the so-called anticubo-octahedron^l or the truncated hexagonal 



bipyrainid.— The essential point is that, since the above polyhedra have several facets (14 
for the cubo- and anticubo-octahedron^l and 20 for the truncated hexagonal bipyramid^^) , 
one expects the morphology of the adsorbed Pt to be sensibly influenced by the diffusivity 
of Pt adatoms/islands through the edges of the Ru nanoparticles. To assess the possibility 
of diffusion of Pt monomers and dimers, we simulate two edges of a Ru nanoparticle using a 
superstructure described in the next section and calculate diffusion barriers of Pt monomers 
and dimers on its (0001) and (1101) facets and across its edges. Note that our model for the 
facet geometries coincide with some of those detected in experiment^ii^ and are also proper 
to the truncated hexagonal bipyramid.— 

The rest of this article is organized as follows: Section [IT] presents the computational 
details, Section IIIII contains our results and provides some discussion about Pt islands on 
Ru(OOOl) (subsection III! Al) and of Pt diffusion on the (0001) facet and across the edges of 
our Ru nanostructure (subsection IIII Bl) . Section HVl summarizes our results and conclusions. 

II. COMPUTATIONAL DETAILS 

Periodic supercell calculations have been carried out within the density functional theory 
(DFT) framework,— 122 using the plane wave and pseudopotential methods^A as embodied 
into VASP (Vienna ab initio Simulation package) ^^^ and with ultrasoft pseudopotentials.— 
We have used a kinetic energy cutoff of 400 eV for the wave functions and 700 eV for the 
charge density to obtain convergent results with sufficient computational accuracy of the 
lattice constant of bulk Ru and Pt. Brillouin zones were sampled with either the (4x3x1) 
or the (3x3x1) Monkhorst-Pack k-point meshes,— depending on the size of the supercell, as 
we will see. Since the main uncertainty of DFT comes from the exchange-correlation poten- 
tial, we have used two different approximations for the exchange-correlation functional: the 
Perdew and Wang generalized gradient approximation (GGA)^^ and the Perdew, Burke, and 
Ernzerhof (PBE) modified GGA— , and compared some of the results obtained using these 
two approximations. To achieve force relaxation of the studied structures, the total energy 
of the system and the forces acting on each atom are obtained after each self-consistent 
electronic structure and minimized by the conjugated-gradient algorithm.^^ At equilibrium, 
forces on each atom are required to be below 0.02 eV/ A. The diffusion barriers for monomers 
and dimers on the Ru superstructure are obtained by the direct dragging method. The 3D 



graphics presented in this work were generated by the Xcrysden program.— 

We note that in the calculations of the geometry and formation energy of Pt islands the 
usage of supercells may introduce contributions from island-island interactions, as a result 
of the imposed periodicity of the system. This is particularly true for the largest islands in 
the 3x4 and 4x4 supercells in which edge-atoms of neighboring islands are as close as third 
nearest neighbors (NN) of the Ru(OOOl) surface. A simple way to estimate this spurious 
interaction is to consider the interaction energy, E^"^^, between two Pt atoms adsorbed on 
Ru(OOOl) as their separation varies and defined as -E™* = E{2Pt/Ru) — 2£'(lPt/Ru), where 
£'(2Pt/Ru) is the adsorption energy of 2 Pt atoms on Ru(OOOl) in a 4 x 4 supercell and 
-E'(lPt/Ru) is the adsorption energy of 1 Pt atom on the same surface and supercell. We 
find that E^"'^ = -0.203, +0.074, and -0.018 eV, as the separation between the two Pt atoms 
increases from 1*** to 2"'^, and 3^'*^ NN bond length, respectively. The interaction between 
3'"'^ NN is thus expected to be 10 times smaller than that of 1*** NN. Accordingly, Eform, as 
calculated in this work, can be reliable up to ± 0.02 eV for the largest islands. 

A. 2D Pt islands on Ru(OOOl) 

The facets of Ru nanoparticles are first modeled by a 5 layer Ru(OOOl) slab, which is 
the surface known to have the lowest energy.— Pt adatoms are placed on only one side of 
the slab. To avoid the interaction between surfaces and Pt adatoms of neighboring periodic 
supercells we have imposed a 15 A vacuum layer, whereas to reduce the interaction between 
deposited Pt islands, the (0001) surface unit cells is extended to either (3 x 4) or (4 x 4) 
structures depending on the island size. The (3 x 4) and (4 x 4) supercells contained 60 and 
80 Ru atoms, respectively, plus Pt atoms forming the island. Their corresponding surface 
Brillouin zone is sampled with a (4 x 4 x 1) and (3x4x1) k-point mesh, respectively. 

B. Monomer and dimer on faceted superstructure 

To model Pt diffusion through the Ru nanoparticles edges we have taken into consid- 
eration edges formed by facets of (0001) and (1101) geometry, which are among the most 
stable Ru surfaces.-*^ We consider a periodic 3D superstructure containing 116 Ru atoms 
and made of a 4-atom wide Ru(OOOl) facet and two Ru(llOl) facets (see Fig. [1]). The con- 



struction of this Ru supercell, which has 7x4 in-plane periodicity, is achieved by stacking 
five Ru(OOOl) layers: two of 7 x 4, one of 6 x 4, one of 5 x 4, and one of 4 x 4 atoms. 
The so obtained edges, on each side of the Ru(OOOl) have different local geometry which for 
convenience are labeled as A and B. Atoms forming edge A (edge B) are contiguous to hep 
(fee) hollow sites of the (0001) facet. The bottom two layers (see Fig. [T]) were not allowed 
to relax to guarantee the stability of the superstructure. We impose a 15 A vacuum layer 
between periodic superstructures along the direction perpendicular to the surface, as in the 
system described previously. The Brillouin zone is sampled with a (2 x 3 x 1) k-point mesh. 
The adsorption energy and diffusion barriers of Pt monomers and dimers are calculated on 
the (0001) and the (TlOl) facets. 

III. RESULTS AND DISCUSSION 

To determine in-plane slab periodicity, we have calculated bulk lattice parameters using 
PW and PBE approximation for the exchange-correlation functional. The bond lengths of 
bulk Ru and Pt are found to be 2.706(2. 744PBE) and 2.77(2.82PBE) A, respectively, while 
the c/a ratio of bulk Ru is found to be 1.585 (PBE) A. 

A. Pt islands on Ru(OOOl) 

We have calculated the optimized geometric structure and energetics of 1 to 5 Pt atom 
islands adsorbed on Ru(OOOl) using the (3 x 4) supercell and that of 1 to 9 Pt atom island 
on Ru(OOOl) using the (4 x 4) supercell. The relaxed structure and total energy of one Pt 
monolayer on Ru(OOOl) has been also obtained. To characterize the energetic stability of a 
given Pt island, we obtain its formation energy, which is defined as: Eform = -E(Ru+Pt) — 
EiYiv) —nE{Piat) ■ Here, £'(Ru+Pt) is the total energy of a Ru slab adsorbed with a n-atom 
Pt island, while -E(Ru) and E{Ftat) denote the total energies of the clean Ru slab and a 
free Pt atom, respectively. Note that the formation energy of stable structures is negative. 
The structure with lowest average formation energy per Pt atom, Eform/n, will thus be 
distinguished as the energetically most favorable one. 

Fig. [2](a) presents Eform/^ as a function of the size of the island, n, for n = 1-5 on the 
(3 X 4) supercell, as provided by PBE and GGA. Fig. [2](b) displays the dependencies of 



Eform/n, for n = 1-4, 6, 7 and 9, on the (4 x 4) supercell. Note that the values n = 12 
and n = IQ represent full monolayer coverage of Pt for the (3 x 4) and (4 x 4) supercells, 
respectively. We find that both PBE and GGA provide similar qualitative dependencies: 
the larger island, the stronger the bonds. Such a trend culminates and is confirmed at full 
monolayer Pt coverage, which provides lowest Ejorm/n. 

The effect of atom detachment (from Pt islands) on Eform has been studied as well. Fig. [3] 
shows two configurations considered for the 7-atom Pt island adsorbed on Ru(OOOl). We find 
that the detachment (transition from the left to the right configuration in the figure) causes 
an increase in Eform from -38.55 eV to -37.52 eV. Similar results have been obtained for 
the islands of other sizes. The increment of energy per Pt atom upon detachment, however, 
does not depend significantly on the island size and vary in the range of 0.11 - 0.14 eV. For 
instance, for a 2-atom Pt island, detachment leads to an increase in Eform from -10.46 eV 
to -10.24 eV while, for a 3-atom island, the energy increases from -15.88 eV to -15.52 eV 
upon the detachment. 

We thus obtain a clear trend: the larger a two-dimensional Pt island (up to 1 ML) is, 
the lower its formation energy per atom. Thus, assuming that the free energy of the system 
in consideration is dominated by its DFT total energy, we conclude that Pt tends to wet 
Ru(OOOl). We have also performed calculations for some 3D Pt islands on Ru(OOOl) and 
found that their 2D isomers have lower energy. For example, two configurations of a 9-atom 
Pt island with 2D and 3D structures (see Fig. Hj) were found to have Eform = -50.06 eV and 
-48.54 eV, respectively. 

The above results indicate that Pt atoms overpower the stress effects derived from the Pt- 
Ru bond length misfit and have a propensity to increase their local coordination on Ru(OOOl) 
regardless of the chemical nature of the bonds, thus forming 2D islands up to ~ 0.56 ML 
and, presumably, up to 1 ML (see Fig. [2]). Of course our calculations do not include the 
case of larger Pt islands (more than 10 atoms) and leave open the question of formation of 
3D structures in such cases. Nevertheless, there is experimental evidence that Pt may form 
2D pseudomorphic islands up to full coverage through vapor deposition.— In turn, from the 
standpoint of ab initio calculations, our conclusion that Pt grows pseudomorphologically is 
justified only at low temperatures and after a long relaxation time. To take into account the 
kinetics of the system, we have looked at the diffusivity of a Pt adatom in the event that it 
be adsorbed on top of Pt islands. Consider, for example, the case presented in Fig. HI The 
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question is thus whether the Pt adatom will be trapped on top of the Pt island (as shown 
in the right inside of Fig. H]) or will it descend at the kink site on the Ru(OOOl) surface 
(left inside of Fig. Hj). Our calculated activation barrier of the Pt monomer diffusion on the 
8-atom Pt island from the hep to fee sites is ~ 0.23 eV and from the fee to hep is ~ 0.09 eV. 
The activation barrier for the step descent for a Pt monomer from the fee to the kink site on 
the Ru(OOOl) is ~ 0.39 eV (and ~ 1.92 eV for the inverse process). Therefore, the Schwoebel 
barrier of this heteroepitaxial step descent of Pt monomers is ~ +0.30 eV, pointing to 3D 
clustering growth. These opposing results hint why the growth mode of Pt on Ru(OOOl) 
sensibly depends on experimental conditions.-"^ Spontaneous deposition of Pt, for example, 
has led to the conclusion that 3D structures follow the 2D pseudomorphic islands if the 
Ru nanoparticles are immersed for sufficiently long periods of time (not specified) into a 
O.IM HCIO4 solution containing O.lniM H2PtCl6,— whereas spontaneous deposition by two- 
minute immersion into 0.1 mM [PtClg]^— and subsequent rinsing with 0.1 M H2SO4 and 
ultra pure water lead to Pt nanoparticles with columnar shape of 10-15 MLs and 3-5 nm.— 
For the purpose of our analysis, however, we shall notice that both vapor deposition below 
0.6 ML4^ and spontaneous deposition via short immersion into Pt-containing solutions^ 
find that the probability for Pt adatom to arrive on top of a Pt island is low. Indeed, the 
Ru nanoparticles that we are trying to model are decorated with Pt at very low coverage 
(~ 0.1 ML) via spontaneous deposition by immersion in a [PtCle]^" and rinsed with 0.1 M 
H2SO4 solution at room temperature,- therefore, only formation of 2D islands is expected. 
Besides, there is indication that, even in the case of positive Schwoebel barriers, as long as 
the atom deposition rate is low and the diffusion speed fast, the reflective property of steps 
is valid only at temperatures well below room temperature for fcc(lll) metals.— 

There still remains the question: why, despite the misflt, 2D conflgurations are favorable? 
In order to grasp further understanding of the issue, one must realize that a number of 
counter-trends must be taken into account. Firstly, the 2D growth is not all that surprising 
for the reason that the bond- length misflt between Pt and Ru amounts only to a few percent. 
Secondly, in the 2D Pt islands under consideration, the number of NN fluctuates from 3 
(single atom) to 9 (full monolayer). Decrease in the number of NN usually causes reduction of 
the equilibrium interatomic distances. Indeed, we flnd that for a free standing Pt monolayer, 
in which every Pt atom has only 6 nearest neighbors, the equilibrium Pt-Pt NN distance 
is much shorter (~2.6 A) than that in bulk Pt (~2.8 A) and even shorter than the Ru-Ru 



NN distance in bulk Ru (~2.71 A). The misfit in low dimensional structures is thus not a 
well defined quantity because of the dependency on the coordination number of the atoms 
in question. Correspondingly, the issue of stress induced by the bond-length misfit between 
the Pt nanostructures and the Ru surface atoms, neither one of which is expected to be at 
the bulk value given the diversity of their local geometric environment, is thus not clear for 
these heteroepitaxial system. 

In addition, Barcaro et al. have pointed out that heteroepitaxial Volmer- Weber (3D) 
growth is favored for those systems in which the adsorbed metal interacts more weakly with 
the substrate than with atoms of the same species.— In our particular case, however, we 
obtain that the Ru-Pt bond may in fact be stronger than that of Pt-Pt since our calculations 
show that the average formation energy per atom of a Pt monolayer on Ru(OOOl) (~ 5.89 
eV) is higher than both the cohesive energy of Pt bulk (~ 5.85 eV) and the average cohesive 
energy per atom of the Pt surface (~ 5.6 eV). 

For the Pt atoms on the top of a hep metal such as Ru, there is also an incommensurability 
in bulk structure. We find that the bulk NN bond length of Pt atoms certainly decreases 
when they arrange in an hep structure. In that case, the bulk bond length misfit between 
Pt and Ru decreases from ~ 2.8 to ~ 1.4 %. Furthermore, the surface interlayer distance 
in Pt(lll) expands to 2.49 A (~ 1.0 % with respect to bulk), while that of the hypothetical 
Pt(OOOl) contracts to 2.39 A notwithstanding that intralayer NN distances are 1.8 % smaller 
than in the fee bulk. Therefore, both the reduced coordination in adsorbed Pt clusters and 
the Pt-Ru bond strength may quench the anticipated effect of bulk bond-length misfit in 
such a way that Pt atoms adsorbed on Ru(OOOl) surface tend to form 2D islands as large as 
possible up to one monolayer. 

B. Modeling Pt diffusion on the (0001) facets and through the edges of Ru nanos- 
tructures 

From the previous section, we have gained an understanding of the tendency of Pt atoms 
to form 2D islands, wetting the Ru(OOOl) surface rather than clustering in multiple 2D or 
3D structures. The experimental evidence, however, suggests that Pt islands maintain small 
sizes on Ru nanoparticles,^ say ~0.5 nm (5 to 7 atoms), for 0.1 ML coverage. The difference 
in the characteristic of Pt surface alloys on Ru(OOOl) and on Ru nanoparticle is part of 
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the reason for the substantial reactivity of Pt adatoms on Ru nanoclusters and not on Ru 
surfaces. If Pt atoms were to make as many bonds as possible on the Ru nanoparticle as 
they do on the Ru surface, one would expect that, even for low coverage (~0.1 ML), a 
large island should totally cover one of the facets of the Ru nanoparticles. For example, a 
hep Ru nanoparticle of 2.5 nm with the proposed anticubo-octahedral structure^i^I could 
hold roughly ~7 (~4) Pt atoms per squared (triangular) facet for homogeneous coverage 
of 0.1 ML (~ 70 Pt atoms), whereas the same coverage coalesced into a single island could 
totally cover one of the squared facets, which seems not to be the case.- Note that a similar 
argument would follow for a truncated hexagonal bipyramid. 

One main difference between the infinite surface and the nanoparticle is that the latter 
exhibits edges dividing its facets. It is natural to assume that they prevent Pt coalescence 
on the Ru nanoparticles. If this is true, 2D islands are formed on each facet, but they do 
not join together into a large unique island because the edges prevent those initial small 
islands from diffusing to other facets, thus persisting as few-atom 2D islands. Support for 
the aforesaid reasoning will be attained below by comparing the barrier for Pt atoms to 
diffuse on a (0001) facet with that for Pt to diffuse across the edges towards a (1101) facet. 

1. Pt Monomers 

The Ru nanostructure used for the above purpose is displayed in Fig. [TJ It possesses 3 
hep (hcpl, hcp2, hcp3) and 3 fee (fcc4, fcc5, fcc6) non-equivalent hollow sites on the (0001) 
facet, as shown in Fig. [5t^a). The calculated adsorption energies {Eads) of Pt monomers on 
these sites are listed in Table [H We find that for all hep sites Eads is higher than for any 
fee site. Table [I] shows in addition that Pt monomers preferably sit on sites surrounded 
by 2 edge atoms (fcc4 and hcp3), rather than on those surrounded by only one (fcc6 and 
hcpl) or none (fcc5, hcp2) edge atom. Note also that the adsorption energy of a monomer 
on the (0001) facet is lower than on the infinite surface (5.13 eV). Across edge A, the first 
(1101) available site is a 4-fold hollow site, denoted by "1" in Fig. [5](b), whose adsorption 
energy, 5.66 eV, is substantially higher than that on hep sites of the (0001) facet. Across 
edge B, the first (1101) available site is a 3-fold hollow site, denoted by "2" in Fig. [5](b), 
whose adsorption energy is 4.92 eV. The above results, incidentally, suggest that there may 
be a propensity of Pt to deposit on (1101) facets in the long range. 
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The calculated diffusion barriers, AE, through edge A, appear to be highly asymmetric 
(see Fig. EJ^a)): AE(hcp3 -^ "1") = 0.49 eV and AE("1" -^ hcp3) = 1.10 eV. On the 
other hand, the shortest path for diffusion through the edge B connects the hcpl site on the 
(0001) facet with the three-fold hollow site on the (1101) facet (fcc4 -^ "2"). The barrier for 
diffusion along this path is found to be 0.22 eV, which is just slightly higher than the barriers 
for diffusion on the facet. However, our calculations indicate that the initial state for this 
process (hcpl site) is metastable since its total energy is 80 meV higher than that for the 
monomer adsorbed on the neighboring fcc4 site. The diffusion of the Pt monomer through 
the edge B is thus a two-step process with diffusion from the hcpl to fcc4 site on (0001) and 
then from that to the tree-fold site on (1101). For the above reason, the overall probability 
of diffusion through the edge B is substantially reduced. Even at room temperature, the 
probability to find a Pt adatom on the fcc4 site is ~ 24 times lower than that on the hcpl 
site. Moreover, the hcpl <-> fcc4 barrier is asymmetric (0.20 eV for hcpl -^ fcc4 and 0.12 
eV to diffuse back), which further reduces the diffusion probability. 

In order to estimate the probability of Pt diffusion predicted by these barriers, we have 
calculated roughly the diffusion rate of Pt monomers. The latter, dominated by the expo- 

AE 

nential of the energy barrier, is given by R = Dqc ''b^ , where ks is the Boltzmann constant, 
T is the temperature, and Dq is the diffusion prefactor— whose typical values are of the order 
of ~ lO^^s"^. By noting from above that the diffusion rate through the edge B is reduced 
by a factor of ~ 24 with respect to that from fcc4 -^ "2", we obtain that the former is 
effectively ~ 6x10^ s^l, thus significantly smaller than that on the (0001) facet (~ 4x10^ 
and ~lxlO^'^s^^, for hcpl -^ fcc4 and fcc4 — ^ hcpl, respectively). On the other hand, the 
diffusion rate through edge A (~ 3 x 10^^ and ~ 6 x 10'^ s" 1, for " 1" — > hcp3 and hcp3 -^ " 1" , 
respectively) is found to be three orders of magnitude lower than that through edge B and 
at least five orders of magnitude lower than that on the (0001) facet. Notice that since the 
barriers to diffuse back from the (1101) facet to the (0001) facet through edge B ("2" — > fcc4) 
is only slightly larger (0.28 eV) than the barrier for diffusion in the opposite direction, diffu- 
sion across edge B may be actually inefficient. The rate at which Pt monomers return to the 
(0001) facet can nevertheless be expected to be lower than that by which they diffuse to the 
neighboring four-fold hollow site on the (1101) facet, since the barriers for the latter process 
are significantly lower. The diffusion from the three-fold hollow site to the four-fold hollow 
site on the (1101) facet is found to be a two-step process. It involves the diffusion to an 
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intermediate local minimuin posing a 0.12 eV barrier and a subsequent step with a barrier 
of only 70 meV. The barriers to return to the three-fold hollow site are, correspondingly, 
0.20 eV and 0.6 eV. These values confirm that, at least energetically, monomers may have a 
propensity to remain on the (1101) facet once they occupy it. However, exact implications 
of these rates and the impact of competing processes can only be visualized after kinetic 
effects are properly included, such as in kinetic Monte Carlo simulations which we leave as 
a further task. In summary, our results indicate that edges compel Pt monomers to remain 
on the facet where they are initially adsorbed, thus preventing the formation of large Pt 
islands on the Ru nanoparticles, although Pt monomers may diffuse across some edges more 
efficiently than across others. Of course, since the work reported here has not considered 
facets of the Ru nanoparticles that exhibit other geometries, the relative adsorption energies 
and diffusion barriers of these may expand and tune the occupancy landscape from the one 
we have inferred. 

2. Pt dimers 

The results of previous sections indicate that clustering of 2D Pt islands on the Ru 
facets occurs readily and that, although the Ru edges hinder the inter-facet diffusivity of Pt 
monomers, they do not unquestionably prevent it at room temperature. Yet, in order for 
Pt islands to coalesce and form larger islands, they would have to diffuse through the edges. 
We therefore turn to the calculation of the energy barrier for a dimer to diffuse from hep to 
fee sites (on the (0001) facet) and through edge B, the easy edge for monomers to cross. 

Some of the sites that Pt dimers may adopt on the (0001) and the (1101) facets, as well 
as the corresponding average formation energies per atom, are shown in Fig. [7l We find 
that Eform/n of the dimer increases (by ~ 0.12 eV) with respect to that of the monomer, 
suggesting that dimers would preferably form rather than diffuse as monomers through the 
easy edges. As in the case of monomers, dimers prefer to sit on hep sites on the (0001) facet 
(see Fig. El^a) and (b)); similarly, the adsorption energy of Pt dimers at hep sites near edge 
B (only one edge-atom neighbor) and in the middle of the Ru stripe is almost the same (see 
Fig. [7](a) and (c)). On the (0001) facet, when one of the atoms in the dimer comes closer to 
the edge and its coordination is reduced from 5 to 4 (compare Fig. [Tl^c) and (d)), E form In 
drops ~ 0.16 eV, suggesting that there is a higher barrier for Pt dimers to approach the 
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edges to the point where its atoms become more undercoordinated. On the (1101) facet, 
analogously, Ejorm/n is 0.52 eV lower for a dimer across the (0001) and (1101) edges (see 
Fig. [7](e)) than for a dimer sitting on the (1101) facet close to the edge (see Fig. Cl^f)). 

As shown in Fig. [HI the barrier for the dimer to diffuse from fee to hep sites has a height 
comparable to that for monomers while it is 2.5 times smaller for the inverse process, from 
hep to fee. The diffusion rates (1.4x10^ and 8.8x10^ s~^, correspondingly) thus differ by 
four orders of magnitude, indicating that diniers are most of the time at hep sites. This 
presents another indication that Pt dimers (or larger islands) would not leave the (0001) 
facet since, for it to diffuse across edge B, it must be on fee sites. 

As shown in Fig. [9], the diffusion across edge B comprises two stages. The initial state of 
the first stage corresponds to the configuration shown in Fig. [7|(d) in which one atom is on 
a hcpl site and the other is slightly beyond the edge whereas, in the final state, the atom at 
the hcpl site diffuses to a fcc4 site and the other moves to a type "2 " site (see Fig. El^b)) of 
the (1101) facet, as shown in Fig. [7](e). The energy barriers for the dimer to diffuse back and 
forth from the initial and final states are shown in Fig. [9] while the corresponding diffusion 
rates for these processes are low (~ 10^). Indeed, the energy barriers of the first stage are 
very similar to those for dimer diffusion from hep to fee sites (see Fig. [8]) and for monomer 
diffusion through edge A (see Fig. [6]). Notice also that the final state described above is only 
an intermediate stage of the diffusion towards the (1101) facet. For the second stage, the 
initial state is naturally the configuration shown in Fig. Wi^) while the final state corresponds 
to that shown in Fig. W(i), in which both atoms sit on the (1101) facet. The energy barrier 
to move from the initial to the final state (see Fig. [9]) produces also a low diffusion rate of 
6x10'^ s~^, while the inverse process, whose barrier is 3 times larger, provides a diffusion 
rate of ~1.4xl0~^^s~^, indicating that dimers like monomers on the (1101) facets will most 
likely remain there. In short, edge B, which offers a low-barrier diffusivity path to monomers, 
renders to Pt dimmers two energy barriers for facet-to-facet diffusion, each representing a 
diffusion rate at least three orders of magnitude lower than those of the monomer. We 
expect the diffusion rates of trimers and other n-mers of Pt on Ru nanoparticle facets to be 
even lower than that found here for dimers. 
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IV. SUMMARY 

We have calculated from first principles the energetics and geometry of Pt islands de- 
posited on Ru(OOOl), as well as the energy barriers for diffusion of Pt monomers and dimers 
through the edges intersecting the (0001) and (1101) facets of a superstructure modeling 
a Ru nanoparticle. We find that the low coordination of Pt atoms composing the islands, 
the strong Pt-Ru interaction, and the hep structure of the substrate promote formation 
of increasingly large Pt islands on Ru(OOOl), possibly up to IML, and avoid the 2D/3D 
crossover. On the other hand, the scenario is quite different when Pt atoms are deposited 
on Ru nanoparticles. In a simple model for (0001) factes with edges connecting to (1101) 
orientations, we concur with experimental predictions that Pt atoms arrange homogeneously 
over the facets of Ru nanoparticles by spontaneous deposition and form 2D islands. We also 
predict that these islands do not coalesce into a large unique island because the edges of 
the Ru nanoparticles prevent monomers and dimers from diffusing to other facets. Our 
calculated barriers indicate that there may be some edges in the Ru nanoparticles for which 
the diffusion rate across-edge is several orders of magnitude lower than the diffusion rate 
on-facet, even for monomers. For those edges that may offer relatively low diffusion barriers 
to monomers, our calculated barriers for dimers suggest that dimers or larger islands, whose 
formation is more probable than the monomer diffusion, remain in the facet where they were 
formed since the diffusion rate across-edge is several orders of magnitude lower than that of 
monomers. 
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TABLE I: Adsorption energy (in eV) of a Pt monomer on various sites of the (0001) and (1101) 
facets of our Ru nanoparticle model (see numbering notation in Fig. [5]) 



0001 hep 




0001 fee 




(1101) 




Label 
(Fig. ED 
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(eV) 
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(eV) 
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(eV) 


1 


4.94 
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4.86 


1 


5.66 


2 


4.93 


5 


4.77 
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4.92 
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5.05 


6 


4.78 
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FIG. 1. (color online) Model of the edges of a faceted Ru nanoparticle exposing a (0001) 
facet and two (1101) facets. Different colors distinguish the five layers parallel to the (0001) 
surface constituting the structure. 

FIG. 2. Average formation energy per atom, Eform/n-, as a function of the size, ra, of the 
island for n = 1 - 5 on the (3x4) supercell (upper panel) and for n =1 - 4, 6, 7, and 9 on the 
(4x4) supercell (lower panel). 

FIG. 3. (color online) Two configurations of a 7-atom Pt island (blue) on Ru(OOOl) 
(grey) showing the detachment of one Pt atom. The configuration in the right panel has 
lower Eform/n than the configuration in the left panel by ~0.14 eV (see text). 

FIG. 4. (color online) Two configurations of a 9-atom Pt island (blue) with 2D(left) and 
3D (right) structure on Ru(OOOl) (grey) used in calculations. 

FIG. 5. (color online) Adsorption sites of Pt monomers (red) on the (0001) facet of the 
Ru nanoparticle model (blue). Numbers "1", "2" and "3" indicate hep sites and "4", "5", 
"6" indicate fee sites, (b) Adsorption sites of Pt monomers on the (1101) facets of the Ru 
nanoparticle model. Numbers "1" and "2" indicate four- fold and three-fold hollow sites, 
respectively. 

FIG. 6. Energy barrier for the diffusion of a Pt monomer (a) across the edge A (see 

Fig. [5]): hcp3 < ^ type "1" site, (b) on the (0001) facet: hcpl < > fcc4, and (c) across the 

edge B: fcc4 < — > type "2" site. 

FIG. 7. (color online) Adsorption sites and Eform/n {n = 2) of Pt dimers (red) on the 
facets of the Ru nanoparticle model (blue), (a) Top view of the (0001) facet showing a dimer 
at hcp2 sites, (b) Top view of the (0001) facet showing a dimer at fcc5 sites, (c) Top view 
of the (0001) facet showing a dimer at hcpl sites, (d) Top view of the (0001) facet showing 
a dimer with one atom at a hcpl site and the other beyond the edge of the same facet. The 
coordination number of the atoms is 5 and 4, respectively, (e) Top view of the (1101) facet 
showing a dimer with one atom at a fcc4 site and the other at a type "2" hollow site, (f) 
Top view of the (1101) facet showing a dimer at a type "2" and type "1" hollow sites of the 
same facet. 

FIG. 8. (color online) The upper three panels from left to right illustrate initial, transition, 
and final states, respectively, of the diffusion of a dimmer (red) from hep to fee sites, as seen 
from the top of the (0001) facet (blue). The lower panel shows the corresponding barriers 
of this process and the inverse, from fee to hep sites. 
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FIG. 9. (color online) The upper five panels from left to right illustrate the two-step 
diffusion of the dimmer (red) across the edge intersecting the (0001) and the (1101) facets 
(blue). First, third, and fifth upper panels are local minimum energy configurations of the 
dimmer and the second and forth upper panels are transition states. The lower panel shows 
the barrier for the dimmer to diffuse back and forth from either local minimum energy 
configuration (see text). 
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